Anisotropic solitons in dipolar Bose-Einstein Condensates 
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Starting with a Gaussian variational ansatz, we predict anisotropic bright solitons in quasi-2D 
Bose-Einstein condensates consisting of atoms with dipole moments polarized perpendicular to the 
confinement direction. Unlike isotropic solitons predicted for the moments aligned with the confine- 
ment axis [Phys. Rev. Lett. 95, 200404 (2005)], no sign reversal of the dipole-dipole interaction is 
necessary to support the solitons. Direct 3D simulations confirm their stability. 
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The creation of two- and three-dimensional (2D and 
3D) solitons is a challenge to experiment in optical and 
atomic physics Both "light bullets" in nonlinear 

optics [2J and 2D/3D matter- wave solitons are objects 
of profound significance. In optics, the best result was 
the creation of partially localized (quasi-2D) spatiotem- 
oral solitons in crystals with quadratic nonlincarity 
As concerns matter waves, quasi- ID bright solitons 
were demonstrated in self-attractive Bose-Einstein con- 
densates (BECs) of 7 Li and 85 Rb atoms and gap 
solitons were created in a self-repulsive 87 Rb condensate 
trapped in an optical lattice (OL) A fundamental 
impediment to the creation of 2D and 3D bright solitons 
with attractive cubic nonlinearity is their instability to 
collapse. Theoretically considered stabilization schemes 
include the use of OLs @ and periodic time modulation 
of the nonlincarity by means of a Feshbach resonance 
(FR) in ac magnetic field [?J • Another approach assumes 
using 2D or 3D OLs to create gap solitons in respective 
settings with repulsive nonlinearity Q. So far, however, 
no experimental realization of matter-wave solitons in a 
2D or 3D setting has been reported, and other stabiliza- 
tion mechanisms are highly desirable. 



One such mechanism may utilize the nonlocal interac- 
tion between dipolar atoms in the quantum Bose gas Q . 
Potential implementations include BECs of magnetically 
polarized Cr atoms permanent electric dipoles in 

quantum gases of heteronuclear molecules 1 11], and elec- 
tric dipoles induced by laser illumination [12j or by strong 
dc electric field. In this context, it is natural to consider a 
nearly 2D "pancake" geometry, the simplest (isotropic) 
configuration having dipole moments aligned along the 
tight-confinement direction, denoted here as z. This ap- 
proach to the creation of quasi-2D solitons in dipolar 
BECs was elaborated in Ref. [l3j, using a variational 
approximation based on the Gaussian ansatz -0i SO = 
7r -3/4 Q ,i/2 7 i/4 exp |_ t a t x 2 + y 2^ + 1Z 2^ / 2 ] ^ w h ere a 

and 7 are the inverse squares of the radial and axial 
widths, respectively. This trial- function generates the 



following Gross-Pitaevskii (CP) energy functional, 
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Here and below, length, time, and energy are scaled as 
r — > r/l±, t — > dJ±t, and E — > E/ (hu>±), where oj± is the 
trap frequency in the z direction, and l± = yJU/mu^. 
The interaction strengths are g = 4itNa s /l± and gd = 

NcPm (h 2 l±) , where a s > is the s-wave scattering 
length, d and m the atomic dipole moment and mass, iV 
the number of atoms, and the wave function is subject 
to normalization J \ip (r)| 2 dr = 1. 

For fixed 7 in the absence of dipole-dipole (DD) inter- 
actions (gd — 0) , both kinetic and contact-interaction en- 
ergies (the first and fourth terms in Eq. (PJ) scale as the 
inverse square of the radial size in the 2D plane, resulting 
in the well-known instability of the Townes solitons [l^ ] . 
Stabilization is shown by fixing 7 = 1 in Eq. ([T]) and 
searching for a minimum of the energy as a function of 
a, which results in the necessary condition [l3j ]. 

(2/3V2V) g d < 1 + . 9 /(2vr) 3 / 2 < - (4/3V2^) g d , (3) 

that may only hold for gd < 0. The interpretation of 
this condition reveals the origin of self-trapping in the 
proposed configuration. Since the DD interaction is re- 
pulsive and attractive, respectively, for "side-by-side" 
and "head-to-tail" dipole pairs, the overall dipolar en- 
ergy varies from negative for a prolate structure, with 
L p ~ a -1 / 2 <C L z ~ 7 -1 / 2 , through zero for the isotropic 
case, L p — L z , to positive for the oblate condensate with 
Lp 3> L z (this variation is evident in the sign change of 
function ([2]) at k = 1). Fixing L z , a maximum may be 
attained for the total energy as a function of L p , provided 
that the DD term is sufficiently large to reverse the sign 
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FIG. 1: (color online) Gross-Pitaevskii energy functional fTJ, 
as a function of variational parameters a and 7, with g — 
500 and (a) g d /g = -0.10, (b) g d /g = -0.20, (c) g d /g = 
—0.40. Weak dipolar interaction in (a) leads to expansion 
of the condensate, whereas strong attraction in (c) results in 
collapse. For intermediate strengths of the dipole interaction 
in (b), stable quasi-2D solitons are possible. 



of the sum of the kinetic and contact-interaction ener- 
gies, all scaling as 1/L 2 . This unstable maximum may 
be turned into a stable minimum if the sign of the dipolar 
interaction strength is reversed to gd < 0. This goal may 
be attained by using rotating external fields [151 ]. but in 
combination with the necessity to reduce contact inter- 
actions in order to make the DD interaction dominant, 
it introduces an essential complication to the experiment 
[l3j |. and it is not possible for electric moments induced 
by the polarizing dc field. 

Condition ([3]) does not guarantee a true minimum of 
E{a, 7) in the (a, 7) plane. For weak DD interactions, 
Eq. © does not hold, and all fixed-7 energy curves are 
monotonic in a, as shown in Fig. [TJa). For larger val- 
ues of |<7d|j a shallow local minimum (elliptic point) does 
appear, as in Fig. [Ijb). It gives rise to a quasi-2D soli- 
ton, which is isolated from collapsing in the direction of 
a, 7 — > 00 by a saddle point. However, with further in- 
crease of the strength of the DD interactions, the elliptic 
and saddle points merge and disappear, resulting in an 
essentially 3D collapse, as seen in Fig. HJc). This is the 
source of the "breathing instability" in Ref. [131 ]. 

Related works considered the stability of vortex lat- 
tices [l6| and Thomas-Fermi-like trapped states [13] in 
dipolar BECs in the same oblate geometry. Also rele- 
vant are theoretically elaborated examples of isotropic 
optical solitons supported by thermal and other nonlocal 
nonlinearities in the 2D geometry (lij ]. 

In this work, we aim to explore a new setting with 
dipole moments polarized in the 2D plane, i.e., ■perpen- 
dicular to the tight-confinement axis, which we now des- 



ignate y, reserving label z for the dipolar axis. With this 
notation, the GP energy functional is 
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Substituting the general anisotropic ansatz, VWso = 
tt" 3 / 4 (a^ 7 ) 1/4 exp [-(1/2) (ax 2 + (3y 2 + 7 z 2 )] , into ©, 
we find (with k x = -yAy/a, n y = \f^f[0) 
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While h(n x , K y ) may be expressed in terms of elliptic 
integrals, its behavior is more readily seen directly from 
Eq. ©. We first consider the variation of the characteris- 
tic width, L z = 1/^/7, along polarization axis z. At small 
L z (7 ct,(3), we have h(K x ,K y ) — * /i(oo,oo) = 2. On 
the other hand, h(n Xl n y ) — > h(0, 0) = —1 for 7 <§; a, (3. 
Thus for gd > 0, i.e., for the natural sign of the DD inter- 
action, Eq. ^ demonstrates a switch from repulsion at 
small L z to attraction at large L z , thus predicting a sta- 
ble bound state as concerns the variation of L z , provided 
that gd is large enough. 

While the above consideration makes stable self- 
trapping along the polarization axis evident, behavior 
along the other unconfined direction, x, is more sub- 
tle. Indeed, one might expect expansion of the con- 
densate in this direction, as the side-by-side DD inter- 
action is repulsive, However, examination of expression 
yields the asymptotic limits, h(0,n y ) = —1 and 
/i(oo, K y ) — (2k v — 1) / (n y + 1), so that for gd > and 
K y < 1/2, the DD interaction remains effectively attrac- 
tive even for arbitrarily large values of L x . For K y > 1/2 
(i.e. for L z > 2L y ), the DD interaction changes from 
the attraction for k x — > to repulsion at k x — > 00. Since 
the attractive part of the DD-interaction energy scales as 
\[u, whereas the kinetic energy in the x-direction scales 
as a, a local minimum can be provided by the interplay 
between them. 

The DD interaction provides for effective attraction at 
large L Xl L z and prevents expansion if it asymptotically 
overcomes the contact interaction. For gd > 0, we thus 
obtain from ([5]) a threshold condition for the existence of 
the soliton, g/gd < 47r/3 « 4.19. In Fig[2l we plot equal- 
energy surfaces and their cross sections in the (a, (3, 7) 
space, for g/gd — 0.911. The expected local energy min- 
imum exists, corresponding to an elongated soliton with 
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FIG. 2: (color online) (a) Equal-energy surfaces around the 
local minimum of Gross-Pitaevskii energy functional in 
the variational parameter space of inverse squared widths 
a,/3, 7, for g = 10 and g/gd = 0.911. Energy values range 
from E = 0.4907 for the innermost surface to 0.4987 for 
the outermost one, in steps of AE = 2 x 10~ 3 . Panels 
(b),(c),(d) depict cross sections through the minimum-energy 
point, (a,/3, 7 ) so i = (0.136,1.08,0.037), in the (0,7), (/3, 7 ) 
and (a, (5) planes, respectively. The color bar in (b) applies 
also to (c) and (d). 
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FIG. 3: (color online) Numerically simulated three- 
dimensional evolution: (a) Density of numerically-exact soli- 
ton solution at y = 0, obtained by the imaginary-time prop- 
agation of the Gaussian profile i/wiso with minimum-energy 
parameters a — 0.136, j3 = 1.08, 7 = 0.037 borrowed from 
Fig. [2] The 2D density profile of this Gaussian function at 
y — is shown in (e). Results of the propagation in real time 
are presented in (b)-(d) and (f)-(h) for the initial functions 
(a) and (e), respectively. Density profiles along axes x,y, and 
z, with perpendicular planes going through the soliton's cen- 
ter, (xo,yo,zo) = (0,0,0), are shown as functions of time. 
Coefficients g and gd are the same as in Fig. [2] 



The existence and stability of the solitons predicted 
by the above analysis have been verified by simulations 
of the underlying 3D GP equation corresponding to en- 
ergy functional ([!]). The soliton shown in Fig. a) was 
obtained by the 3D propagation in imaginary time of 
the Gaussian minimum-energy ansatz depicted in Fig. 
[3te). Forward propagation in real time, starting from 
this soliton solution (Figs. EUb)-(d)), demonstrates per- 
fect self-focusing. The stability is also demonstrated by 
the real-time 3D propagation starting with the Gaussian 
approximation, as shown in Figs. [2If)"( n )- The observed 
long-lived, slowly damped oscillations indicate the exci- 
tation of an internal mode [20l | , which is typical to stable 
solitons in nonintegrable systems. 

An interesting issue, which will be considered sepa- 
rately, is interactions between the solitons. Collisions 
between isotropic solitons supported by the inverted DD 
interactions were simulated in Ref. [ijj], demonstrating 
their fusion into a single soliton. In the present setting, 
anisotropy is expected to produce various outcomes of 
collisions, depending on the angle between the soliton's 
polarization and collision direction. 

Another dynamical behavior, unique to anisotropic 
solitons, may be attained by rotation of the external mag- 
netic field used to align the dipole moments. Whereas the 
slow rotation of the field is adiabatically followed by the 



soliton, as shown in Fig. HJa), it lags behind faster re- 
volving fields, see Fig. 0Jb). For instance, at u±t = 10 
the field is directed along x, whereas the long axis of 
the density profile is not parallel to it. This drag is ac- 
companied by significant deformation and excitation, as 
evident from the expansion of the density distribution in 
Fig. 0Jb). The amount of excitation around the rotating 
soliton (deviation from adiabaticity) smoothly increases 
with the rotation velocity, roughly like in the case of Lan- 
dau - Zener tunneling. For slow rotation, it is is clearly 
seen to be exponentially small in the rotation period. 

To estimate the experimental feasibility of the pro- 
posed quasi-2D solitons in dipolar BECs, we notice that 
the above-mentioned threshold for the existence of the 
anisotropic solitons is gd/g > 0.24. For 52 Cr with 
magnetic dipole moment d — 6/is and s-wave scatter- 
ing length a « 100ao {fiB and ao are the Bohr mag- 
neton and radius, respectively), this ratio is gd/g — 
fi^cF 'm 1 '(167T 2 ?! 5 a) = 0.036 [l0( (/j,o is the vacuum perme- 
ability). Therefore, attenuation of the direct interaction 
between atoms may be necessary. A recent experiment 
has demonstrated that this interaction may be practically 
switched off by means of the FR in the condensate of Cr 
atoms, making the DD interactions absolutely dominant 
2l| . Alternatively, using BECs of dipolar molecules, 



with the electric dipole moment on the order of Debye, 
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FIG. 4: (color online) The rotation of the soliton, driven by 
the polarization field revolving at angular velocity f2, start- 
ing with the anisotropic soliton shown in Fig. EJa), for 
Q = 7rjx/200 (a) and = nu±/20 (b). The interaction 
strengths are the same as in Figs. [2] and [3] 



will certainly provide sufficiently strong DD interactions. 

In conclusion, we have demonstrated, using the varia- 
tional analysis and direct simulations of the GP equation 
in 3D, that the interplay of the ordinary repulsive contact 
interactions and long-range DD forces may give rise to 
stable quasi-2D anisotropic solitons in the dipolar BEC, 
if atomic moments are polarized in the 2D plane. Unlike 
the previously considered isotropic solitons polarized per- 
pendicular to the plane, the stability of the anisotropic 
solitons does not require artificial inversion of the sign 
of the dipole-dipole interaction. The anisotropic soliton 
may adiabatically follow slow rotation of the polarizing 
field. 
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